Resonant acousto-optics in the terahertz range: TO-phonon polaritons driven by an 

ultrasonic wave 
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The resonant acousto-optic effect is studied both analytically and numerically in the terahertz 
range where the transverse-optical (TO) phonons play the role of a mediator which strongly couples 
the ultrasound and light fields. A propagating acoustic wave interacts with the TO phonons via 
anharmonic channels and opens band gaps in the TO-phonon polariton energy dispersion that 
results in pronounced Bragg scattering and reflection of the incoming light. The separation in 
frequency of different Bragg replicas, which is at the heart of acousto-optics, allows us to study 
the resonant acousto-optic effect in the most simple and efficient geometry of collinear propagation 
of electromagnetic and ultrasonic waves. The acoustically induced energy gaps, Bragg reflection 
spectra, and the spatial distribution of the electric field and polarization are calculated for CuCl 
parameters, in a wide range of frequencies and intensities of the pumping acoustic wave. Our results 
show drastic changes in terahertz spectra of semiconductor crystals that opens the way for efficient 
and accessible manipulation of their infrared properties, by tuning the parameters of the acoustic 
wave. 

PACS numbers: 71.36.+C, 43.35.Sx, 42.70.Qs 



I. INTRODUCTION 

Acousto-optics is a well established discipline and has 
a relatively long history which started in 1922 from a 
pioneering work by Leon Brillouin, 1 who predicted light 
diffraction by an acoustic wave (AW) propagating in a 
medium. Ten years later this prediction received clear 
experimental confirmation,^ however, the dielectric con- 
stant modulation caused by the photoelastic effect turned 
out to be extremely small (of order of 10~ 5 ). Such a 
modulation leads to a similarly small light diffraction ef- 
fect, which is normally seen only in the lowest diffrac- 
tion order. This picture is typical for conventional, non- 
resonant acousto-optics^ In contrast, resonant acousto- 
optics, which has been introduced very recently^— deals 
with acoustically induced spectral changes which are 
comparable to the background signal. This is achieved by 
a resonant coupling of the incoming light to such funda- 
mental excitations in the solid state as excitons or trans- 
verse optical (TO) phonons leading to the formation of 
polaritons^ These excitations, in turn, are strongly cou- 
pled to the AW via the deformation potential (in case 
of excitons) or phonon anharmonicities (in case of TO 
phonons). This results in drastic changes of the opti- 
cal properties of solid state structures when they are 
exposed to an ultrasonic wave. The predicted strong 
resonant acousto-optic effect has been recently observed 
in semiconductor microcavities parametrically driven by 
a surface AWi 10 ' 11 A similar effect has been proposed 
and observed in photonic crystal structures, in which 
the distributed feedback cavity enhances the light-sound 
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A travelling AW produces a propagating periodic per- 
turbation in a solid state structure. To some extent, such 
periodicity is similar to that in a photonic crystal. Thus 



the Bragg diffraction of light typical for photonic crystal 
slabs is expected also in structures with travelling modu- 
lations. In the case of an acoustically induced grating the 
periodic modulation propagates in space with the speed 
of sound. Therefore we are dealing with a dynamical pho- 
tonic crystal structure which is essentially different from 
conventional photonic crystals. Apart from a perfect pe- 
riodicity induced by the AW, an obvious advantage of 
this dynamical structure is the ability to tune the Bragg 
lattice constant by changing the frequency of the pump- 
ing AW. A further advantage stems from a fundamental 
property of these dynamical structures, that the Bragg 
replicas, which can be spatially resolved in the far field, 
also acquire a spectral resolution. In this case the scheme 
for the detection of Bragg scattering can be greatly sim- 
plified by allowing the light and the AW to propagate 
collinearly, therefore, using only the spectral/temporal 
resolution. We employ such a scheme in the present pa- 
per assuming both light and ultrasound waves propagat- 
ing in the same/opposite direction that allows us to re- 
duce the strong-coupling theory of a three-body problem 
to an effective one-dimensional dynamical problem. 

The purpose of this paper is to study in detail the 
far-infrared properties of TO-phonon polaritons driven 
by an ultrasonic wave, by solving (numerically) exactly 
a pair of Maxwell's and material macroscopic equations. 
A CuCl semi-infinite crystal is taken as an example in all 
simulations presented in the paper. Since the TO-phonon 
resonant energy typically ranges from several to tens of 
milli-electron- volts (20.3 meV in CuCl) we concentrate 
in this paper on polariton properties in the far infrared 
(terahertz) range. 

We show, in particular, that owing to the anharmonic 
coupling to TO phonons, a coherent ultrasonic excitation 
gives rise to travelling periodic modulations in semicon- 
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ductor structures and to dynamical Bragg scattering of 
electromagnetic waves. In unbound semiconductors, we 
describe this scattering in terms of polaritonic partial 
plane waves which include all important Bragg harmon- 
ics and calculate the quasienergy spectrum of such exci- 
tations. The translational invariance of the quasienergy 
spectrum brings into resonance different polariton states 
which strongly couple to each other. As a result of such 
anticrossing of different polariton branches, band gaps 
induced by the coherent AW open in the quasienergy 
polariton spectrum and manifest themselves as maxima 
in the reflection which are in turn efficiently tunable by 
changing parameters of the pump AW. Expanding the 
electromagnetic field and the TO-phonon polarization 
into partial polaritonic waves of unbound system, we cal- 
culate exactly (with a given accuracy) the full reflectiv- 
ity spectra as well as frequency selected particular Bragg 
replicas and their interference and analyze in detail their 
dependence on the AW parameters. 



We have also developed a perturbative analytic ap- 
proach which allows us to calculate acoustically induced 
band gaps and Bragg reflection spectra and analyze them 
in terms of phonon-assisted transitions between differ- 
ent TO-phonon polariton states. Comparing analytic re- 
sults with the full numerical calculation of the reflection, 
we demonstrate that the analytical approximation works 
well for low and moderate excitation power of the AW. 
Based on the translational symmetry of the AW-driven 
polariton quasi-energy spectrum, the analytic approach 
allows us to show explicitly that the acousto-optic effect 
is resonantly enhanced in down-converted Bragg repli- 
cas. The maximum enhancement occurs at the frequen- 
cies satisfying the Bragg condition^ due to the polariton 
scattering by the AW accompanied by emission of acous- 
tic phonons. 



The paper is organized as follows. Section II presents 
the general formalism of the approach, introducing the 
macroscopic equations valid for collinear propagation of 
the light and ultrasound waves and their solutions in 
terms of plane waves. The quasi-energy spectrum of bulk 
TO-phonon polaritons driven by an AW is calculated and 
analyzed in Sec. III. The vacuum/semiconductor inter- 
face is introduced in Sec. IV where the frequency-resolved 
Bragg replicas and the full reflectivity spectrum are cal- 
culated, both analytically and numerically. In Sec. V, we 
discuss in more detail how the terahertz properties of 
semiconductors change with acoustic frequency, paying 
attention to the limits of very high and very low acoustic 
frequencies as well as to decay and redistribution of the 
electric field. Results are summarized in Sec. VI. Some 
important details of the analytic approach and numerical 
methods used in the present work are given in Appendices 
A-E. 



II. FORMALISM 

The propagating AW is excited through the piezoelec- 
tric component of the electron-phonon interaction. At 
the same time, it couples to TO-phonons via anharmonic 
channels. The lowest order effect is the cubic phonon 
anharmonicity, which is described by the following in- 
teraction, here written in terms of phonon displacement 
operators : 18 ! 19 

V = J2 * (3) (k, p)(ftp + 6t p )(6k-p + &l k+p )(o-k + 4) . 
k P 

(1) 

where and 6k are the acoustic and TO-phonon anni- 
hilation operators, respectively. In case of a propagating 
ultrasonic pump wave, the acoustic-phonon operators in 
Eq. (fT]) are replaced by their classical expectation values 
giving 

V = aJ2b p bl K+p e int + H.c, (2) 
p 

where off-resonant terms are neglected (they are however 
taken into account below in the macroscopic equations). 
The coupling strength 

a cx V4^$ (3) (3) 

is proportional to the AW amplitude (J ac is the acoustic 
intensity) and to the Fourier transform of the third-order 
derivative of the interatomic potential $^ 3 - ) . A careful 
investigation 8 - shows that, to good approximation, any 
dependence on phonon wave vectors in $^ 3 ^ can be safely 
neglected. The operating acoustic frequencies are all lim- 
ited to the linear regime of the acoustic phonon dispersion 
given by 

n(K) = v s K, (4) 

where K is the acoustic wave vector (K = |K|) and v s is 
the sound velocity. 

The macroscopic polarization P(r, i) due to the TO 
phonons is introduced as the inverse Fourier transform 
of the expectation value of the TO-phonon displacement 
operator, (b p (t) + &_ (£)). Then, using the fact that 
D = £fcE + 47rP, where E is the electric field and Eb is 
the background dielectric constant, the TO-phonon po- 
lariton in the field of a coherent pump AW is described by 
a system of coupled Maxwell's and material equations^ 
In what follows, we concentrate on collinear propagation 
of the acoustic and light waves, normal to the semicon- 
ductor surface. In case of a semiconductor occupying half 
space z > 0, the system can be effectively reduced to a 
one-dimensional geometry and is described by the follow- 
ing differential equations (h — 1 is used for simplicity of 
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P(z,t) (5) 



d 2 d 

^ + 27— + u? t + 4crw t cos(Kz - fit) 
ot at 



P[z,t) 



= ^M), (6) 

where uj t and 7 are, respectively, the frequency and 
the damping of bulk dispersionless TO phonons, ujr 
is the Rabi frequency which characterizes the coupling 
strength between light and TO phonons. Equations 
([5]) and © generalize the famous polaritonic Hopfield 
model^ by adding the interaction with the acoustic field. 
They can be derived with the help of Heisenberg equa- 
tions of motion from the Hamiltonian of TO-phonons 
polaritons j 20 ' 21 An essential difference to the conven- 
tional polariton model is the appearance of an oscillating 
term in Eq. ©, which originates from the TO-phonon 
coupling to the coherent AW described by Eq. ([2]) . 

We solve Eqs. ([5]) and © in the semiconductor region 
in terms of partial plane waves. Due to the oscillating 
pump acoustic field, the electric field and polarization 
of each partial wave can be expanded into the following 
Fourier series 



E{z,t) = e^-^J^En 
n 

P[z,t) = e i(kz - u V P n e 



jn(Kz-Qt) 



in(Kz-Qt) 



(7) 
(8) 



where uj and k are, respectively, the leading frequency 
and wave vector of the partial wave of the electro- 
magnetic field. The acoustically-induced perturbation is 
periodic both in time and space and thus it is formally 
required to make a double Fourier transform. However, 
as shown in Appendix [5] this leads to linearly dependent 
combinations of plane waves, so that the double Fourier 
transform of E and P is effectively reduced to the form 
of Eqs. (|7|) and flSJ). We substitute them into Eqs. jSJ and 
(|6]) and equate coefficients at the same exponentials. This 
results in the following matrix eigenvalues problem: 



(k + nK) 2 - J( w + ntl) 2 ] E n = ^r(uj + nQ) 2 P n , (9) 



ui 2 — 2ij(u> + nQ) — (u + nQ)' 



Pn 



2 

+2auJ t (P n+1 + P n _i) = £ -±pt En 



(10) 



We solve the set of coupled linear equations Eqs. ([9]) and 
(j!0[) and find the polariton dispersion ui{k) and the cor- 
responding eigenvectors E n and P n which are the am- 
plitudes of the partial plane waves of the TO-phonon 
polariton propagating in a semiconductor modulated by 
the AW. 

The matrix eigenvalue problem formulated in Eqs. @ 
and (fTUj) is in general complex, and the values of the 



light frequency uj and/or momentum k can be complex 
even if the damping 7 is neglected. Here we distinguish 
two different physical problems: (i) A modified polariton 
quasi-energy spectrum where for any given real value of 
the polariton momentum k we find the complex polari- 
ton energy ui, as done in Sec. lIIII below, (ii) The far in- 
frared properties of the acoustically-driven TO-phonon 
polaritons, such as reflection, scattering, and transmis- 
sion, therefore we are dealing with the real frequency u> 
of the incoming light. In case (ii), considered in Sec. II VI in 
detail, the wave vector k is complex- valued as it describes 
propagation and attenuation of the polariton waves in the 
medium. 

In all numerical simulations done in the present work 
we use the following parameters of bulk CuCl: 22 TO- 
phonon energy w t = 20.28 meV, Rabi splitting ojr = 
14.53 meV, low temperature phonon dampin g 23 ' 24 7 = 
0.2 meV (at T = 5— 10K), background dielectric constant 
Eb = 5.2, and sound velocity v s — 2.02 • 10 3 m/s. The 
coupling constant a, one of the important parameters of 
the model, scales with the acoustic intensity according to 
Eq. ([3]) , with the value of a = 1 meV corresponding^ to 
the acoustic intensity I ac = 14.8kW/cm 2 . 



III. QUASI-ENERGY SPECTRUM OF 
ACOUSTICALLY DRIVEN TO-PHONON 
POLARITONS 

A. Eigenfrequencies 

In unbound semiconductors, the polariton quasi- 
momentum k is a good quantum number which takes 
real values. For nonzero damping 7, the polariton eigen- 
frequencies, in turn, are always complex. Let us however 
concentrate on propagating solutions and disregard, for 
the moment, energetic band gaps where propagation is 
not allowed. Then, assuming 7 = 0, all eigenfrequencies 
also become real. This slightly idealized picture of un- 
damped polaritons in unbound semiconductors will help 
us, however, to understand the basic properties of the 
polariton dispersion modified by the AW. 

Treating k as a real parameter, we solve the eigen- 
value problem, Eqs. (J9)) and (|10[) . finding polariton eigen- 
frequencies uj = io\(k) which correspond to the energy 
bands A = 1,2,... separated by acoustically- induced 
band gaps. Technically, we truncate Eqs. (JSJ and (JTUJ) 
at — n max $J n n max introducing the matrix dimen- 
sion M = 2n max + 1. The system of two Mx M matrix 
equations, containing nonlinear dependence both on fre- 
quency and wave vector is reduced, in Appendix [B] to 
the following linear matrix eigenvalue problem 



WX = wX , 



(11) 



solved by means of complex matrix diagonalization, 
where W is an w-independent 4Mx4M hypermatrix and 
X is a 4M-hypervector containing the Fourier compo- 
nents of the electric field and polarization, defined in 
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FIG. 1: Quasi-energy spectrum of UP (blue) and LP (green) 
branches of undamped TO-phonon bulk polaritons driven by 
an AW with v = Sl/(2n) = 100 MHz (K = 0.3109 ^m" 1 ) and 
a — lmeV. The spectrum without AW, u = uj°(k), is given 
by thin dashed curves. Thick dashed horizontal line shows an 
artificial cut-off in the quasi-continuum of LP states. Panels 
(b) - (e) highlight particular areas in the spectrum. 



spectrum, just in the center and on the borders of the 
Brillouin zone, see Fig. 1 (b), (d), and (e). These gaps 
are acoustically induced spectral regions where polari- 
ton propagation is not allowed and thus they crucially 
affect the terahertz properties of the polaritons, such as 
light reflection and Bragg scattering. The positions of 
the gaps are determined mainly by the frequency fi (or 
wave vector K) of the pumping AW whereas their widths 
strongly depend on both the wave vector K and the in- 
tensity (coupling strength a) of the AW, as discussed in 
detail in the next section. 

In the polariton bands, which are the regions where po- 
lariton propagation is allowed, the modified spectrum de- 
viates from the bare one. This is illustrated in Figs. 1 (d) 
and (e). This deviation become dramatic close to and 
above the TO-phonon resonant energy. Unlike the bare 
polariton dispersion in which the LP branch is bounded 
by the TO-phonon frequency, the AW-induced branches 
of the quasi-energy spectrum, even those which origi- 
nate from the LP branch, can be found at any frequen- 
cies. Thus, they are in principle unbound. They form a 
sort of a quasi-continuum of polariton states detailed in 
Fig. 1 (c), which is due to the acoustically-induced anti- 
crossing and level repulsion. The optical strength of these 
quasi-continuum states, however, drops quickly with en- 
ergy and contributes to the reflectivity only at the fre- 
quencies around the low-energy edge of the Restrahlen 
band. 



B. Acoustically-induced band gaps 



Eqs. (|B8| ) and (|B9J) . respectively. Note that the extra 
factor of 4 in the matrix dimension AM comes from two 
different sources. The existence of two branches in the 
bare polariton spectrum, the so-called upper polariton 
(UP) and lower polariton (LP) branches provide a factor 
of 2. The fact that all eigenfrequencies come in pairs 25 
(oj\ and — wt) provides the remaining factor of 2. 

Clearly, the quasi-energy spectrum ui\(k) has the 
following translational property valid within the same 
branch A: 



£j\{k + sK) — id\(k) + sfl . 



(12) 



where s is an integer. The corresponding eigenfunctions 
are also invariant in a way that, for the same A, 



E n {k + sK) =E n+s {k). 



(13) 



and a similar expression is valid for P n . 

The translationally invariant polariton spectrum &\(k) 
is shown in Fig. 1 as a function of the polariton momen- 
tum k (measured in units of the acoustic wave vector 
K). Due to the spatial and temporal periodicity im- 
posed by the pumping AW, the bare polariton spectrum 
(dashed curves) is folded into the first Brillouin zone and 
then translational shifted, both horizontally and verti- 
cally, into all other Brillouin zones. The coupling of po- 
laritons to the AW opens gaps in the folded polariton 



The band gaps are spectral regions forbidden for po- 
lariton propagation, therefore one can anticipate an en- 
hancement of reflection at the frequencies of incoming 
light close to or inside the gaps. The analytic calcula- 
tion of acoustically induced band gaps, which we show in 
this section, is particularly important for the understand- 
ing of this enhancement and other terahertz properties of 
acoustically driven polaritons, considered later in Sec. lIVI 

The tridiagonal form of the matrix in Eq. (|10j) allows 
us to evaluate the widths of the energy gaps analytically. 
Expressing P n from the diagonal Eq. (|9]) and substituting 
it into Eq. (fT0|) results in 



a n D n E„ - a(a n+ iE n+1 + a n _i.E n _i) = (14) 

with 

D n (k,uj) 



1 

2^t 



uj t — 2vy{u) + nfl) — (u + nil) 



a n (k,u) = — 



a n (k,u>) 
c 2 ( k + nK N 



- 1 . 



(15) 
(16) 



Equation fll4[) can be solved recursively. Introducing 



5 



= E n /E n -i, it is transformed to 



OfnAi - era 

n+lsn+1 



(17) 



With a truncation E-n-1 = E\ = 0, Eq. (fT4|) is solved 
up to the iV-th order. The ratio <r/D n ~ 2a/uj t serves 
as a smallness parameter in our perturbation approach 
which breaks down when the frequency u approaches ei- 
ther the TO-phonon resonance uj t or the bottom of the 
UP branch, ^Juf + wjj. 

The above truncation is equivalent the following 
boundary conditions added to the recursive Equa- 
tion (|T?|): 

= l/£o = 0. (18) 

Taking the first boundary condition £-n-i = as a star- 
ing point, we find £_jv, £-^+1, and so on, ending up 
with the relation ao-Do = cra_i£_i, in which the sec- 
ond boundary condition is used. The above procedure 
turns out to be particularly useful for a perturbative ap- 
proach to the resonant acousto-optic effect discussed in 
Sec. lIVDl where the Bragg replicas are calculated analyt- 
ically in lowest orders. In the present section we evaluate 
perturbatively only the acoustically induced band gaps, 
again assuming 7 = 0. 

In the N = case, Eqs. (fT7|) and (fT8j) result in a single 
equation 

D (k,u)=0, (19) 

which yields the bare polariton dispersion u> — uj (k) con- 
sisting of the LP and UP branches (see dashed curves in 
Fig. 1), with no effect from the AW. The acoustically in- 
duces gaps appear in the polariton spectrum in first and 
higher perturbation orders at specific values of the polari- 
ton wave vector given (neglecting the acoustic frequency 
f2) by the famous Bragg condition 



k N = NK/2, AT =1,2 



(20) 



which corresponds to centers and borders of different 
Brillouin zones (obviously all of them can be translated 
into the first Brillouin zone). The energy positions of 
the gaps, ojjv, are given (in zeroth order) by the bare 
polariton dispersion, 



UJ N 



J°(k N ), 



(21) 



or, in terms of Eq. (|19[) . can be found by solving the equa- 
tion _Do(fcjv ,wjv) = 0. The widths of the gaps, Ajv, are 
evaluated in lowest order in Appendix [D] with the fol- 
lowing result: 
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FIG. 2: (a) Analytical (curves) and numerical (symbols) cal- 
culation of the acoustically induced band gaps in the LP 
branch vs the acoustic frequency v = Q/(2tt) (top axis) and 
the detuning uit— wi (bottom axis), (b) Energy gaps in the UP 
branch vs acoustic frequency (bottom axis) and the detuning 
oji — v/c^f + u% (top axis). In both cases (a) and (b), Ui is 
the energy position of the lowest gap, either in the LP (a) or 
in the UP part (b) of the spectrum. The coupling strength 
a = 0.4 meV, the TO-phonon damping 7 = 0. 



where 



S(u) 



2wwtwfj 



(23) 



(22) 



The gap Ajy opens in the iV-th perturbation order as 
a result of anticrossing of different polariton dispersion 
branches. The anticrossings occur due to the transla- 
tional invariance of the spectrum which brings into res- 
onance states originating from negative-A: and positive-fc 
parts of the bare polariton dispersion, e.g. by a trans- 
lational shift in the (k, w)-plane of the states in the 
negative-fc part of the dispersion by the vector (NK,Nfl), 
i.e. by the same integer of the wave vector and the fre- 
quency of the AW. In other words, the N-th gap opens 
due to the A^-acoustic phonon transition, induced by the 
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pumping AW, which couples two point in the bare disper- 
sion uj (k) having the opposite values of k as given by the 
Bragg condition Eq. (|20|) . As a result, the width of the 

gap A at is proportional to a N (and to I&J 2 ). The form- 
factor S(ui) grows by increasing the TO-phonon compo- 
nent and reducing the light component of the polariton, 
hence the gaps become wider. In fact, since TO phonons 
play the role of a mediator between the light and acoustic 
fields (in our model, these two fields do not interact di- 
rectly) , and thus the more the component of TO-phonons 
in the polariton field, the stronger the resonant acousto- 
optic effect. 

Figure 2 presents a comparison of the analytical and 
numerical calculations of the band gaps at a moderate 
acoustic intensity 7 ac = 2.4kW/cm 2 (er = 0.4meV). The 
analytic results for the widths of the two lowest band 
gaps, Ai and A2 (lines), are in excellent agreement with 
the full numerical calculation (symbols), practically in 
the whole range of the acoustic frequencies and both in 
the UP and LP parts of the spectrum, as seen from Fig. 2. 
At smaller detuning which corresponds to higher (lower) 
acoustic frequency in the LP (UP) branch, the agreement 
gets worse because the analytical model, which is based 
on the assumption that the detuning is large compared 
to the coupling strength, breaks down in this case. At 
the same time, higher-order band gaps like A3 and A4 
are not well reproduced by the analytics. This can be 
explained by the fact that the present analytic approach 
neglects the acoustic frequency Q, assuming that it is 
small compared both to the light frequency and to the 
energy gaps (for details, see Appendix ICj). This, however, 
is not the case of the third- and fourth-order band gap, 
as their widths fall into sub-/ieV-range and thus become 
comparable to or even smaller than f2. 



IV. BRAGG REFLECTION OF 
ELECTRO-MAGNETIC WAVES 

A. Polariton quasi-momenta 

The quasi-energy spectrum and the band gaps dis- 
cussed in detail in Sec. lIIII can be probed in experi- 
ments by a monochromatic electromagnetic wave with 
frequency oj. In this case, the boundary between the 
semiconductor and the vacuum becomes an essential part 
of the system, so that the momentum across the surface 
is no longer conserved. However, for the given frequency 
uj of incoming light, the polariton waves can have only 
specific values of the complex momentum, k = kj(u>). 
They are known in the literature as forced harmonic 
solutions. 26 The complex quasi-momenta Kj(uj) can be 
found as eigenvalues of the same Equations (|9|) and (flO|) . 
In this case we solve the inverse problem to Sec. lIII Al 
the complex quasi-momenta kj(u>) are found as eigen- 
values. The electric field and polarization of the prop- 
agating polariton waves are the corresponding normal- 
ized eigenvectors, E n j and P n j, where an extra index 
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ti(0 (meV) 

FIG. 3: (a) Real part and (b) imaginary part of the po- 
lariton quasi-momentum x(w) calculated with [black and 
blue (dark gray) curves] and without AW [red (light gray) 
curves]. The insets highlight particular regions correspond- 
ing to the lowest acoustically induced gap. The TO-phonon 
damping 7 = 0.2 meV and the AW frequency v — 50 MHz 
(K — 0.1555 fim -1 ). 

j = 0, ±1, ±2, . . . is added in order to distinguish differ- 
ent partial waves. For each frequency, Maxwell's bound- 
ary conditions determine a unique combination of the po- 
lariton partial waves and finally the reflection and trans- 
mission properties which will be discussed in subsequent 
sections. The partial waves themselves are however not 
fully independent. For any integer number s, the wave 
vectors and amplitudes of the partial waves obey the fol- 
lowing translational relations [cf. with Eqs. fH?)) and (|13l) ] 

kj +s (uj) = Kj(u — sf2) + sK , (24) 
E n . j+S (uj) = E n+S j(w - sfi) , (25) 
P n ,j+ S {uj) = P n+Stj (u - sty (26) 

following from the spatio-temporal periodicity of the sys- 
tem pumped by an AW, periodic in time and space. 

The quasi-momenta kj are calculated in the following 
way. The infinite-matrix problem Eqs. ([9]) and (jTDJ) is 
again truncated to — n max ^ n ^ n max [M = 2n max + 1). 
The amplitude P n is expressed from Eq. © and substi- 
tuted into Eq. (flT)]) . The latter equation is then linearized 
with respect to Kj . CIS done in Appendix [B] and brought 
to the form 

VY = kY, (27) 

with a 2M x 2M hypermatrix V and a 2M-hypervector 
Y defined in Eqs. (|Bl9|) and (|B20j) . respectively. There is 
now only one extra factor of 2 in the matrix dimensions 
2M . This is due to the symmetry k — > -k of the ini- 
tial eigenvalue problem. However, there is no more clear 
separation of the eigenvalues into LP and UP branches. 

Technically, we calculate the eigenvalues kj for quite 
big matrices (with M = 161, so that -80 < j < 80), for 
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dashed lines. The real part of the momentum shows the 
typical anticrossing behavior while the imaginary part, 
strictly absent at 7 = a = 0, demonstrates a pronounced 
band gap. Polariton waves for which the frequencies fall 
in the gap region propagate only into a limited area due 
to their finite decay in space, which is determined by the 
imaginary part of the wave vector. With finite damping 
(which is still smaller than the band gap width), the sit- 
uation does not change much, but the imaginary part of 
the wave vector remains finite everywhere in this case, 
see Fig. 4(b). 

B. Bragg replicas: Numerical approach 



FIG. 4: Real and imaginary parts of the polariton quasi- 
momentum h{oj) calculated around the acoustically induced 
gap (a) for zero damping 7 = and (b) for finite damp- 
ing 7 = 0.2 meV, with (solid lines) and without AW (dashed 
lines). Black and red (gray) colors correspond to positivc- 
k unshifted and negative-fc one TA-phonon shifted polari- 
ton dispersion branches, respectively. The AW frequency 
i/ = 50MHz. 

different values of 10 and then select only the central root 

k(lo) = k J=Q (u>) , (28) 

i.e. the root which satisfies the equation x(0) = 0. The 
central root is chosen as best calculated. In fact, it is the 
most distant solution from the boundaries of the trun- 
cated array, so that any distortions of the spectrum due 
to such truncation are minimized in the center of the 
array. We then use this central root >c to improve cal- 
culation of all other wave vectors kj (j 7^ 0) exploiting 
the translational property Eq. (1241) and making a linear 
interpolation Eq. (|B21|) between neighboring frequency 
points. 

Both the real and imaginary parts of the central root x 
are shown in Fig. 3, in the presence of AW and without it. 
Apart from very strong changes in the polariton momen- 
tum slightly below the TO-phonon resonant frequency 
(20.3 meV) and in the whole Restrahlen band above it, 
tiny changes occur near the acoustically induced gaps and 
grow with increasing coupling strength cr, see the insets 
to Fig. 3. These small features, however, are very impor- 
tant for the infrared properties of the system and thus 
we consider them in more detail. 

The real part of the quasi-momentum [Fig. 4(a), top 
panel], if it were calculated for v = 100 MHz, would re- 
produce the quasi-energy spectrum shown in Fig. 1(e) for 
zero damping. The only difference would be the hori- 
zontal connection line in Fig. 4(a). However, the results 
shown in Figs. 3 and 4 are calculated for half that acous- 
tic frequency, which has moved the gap to energies nearly 
half that in Fig. 1 (e) . The energy gap is due to the acous- 
tically induced strong coupling between the two bare po- 
lariton branches crossed in the center of Fig. 4(a), see 



Let us now consider the effects of reflection, trans- 
mission and Bragg scattering due to a semiconduc- 
tor/vacuum boundary at z = 0. In the semiconductor 
area [z > 0), the electric field and TO-phonon polariza- 
tion are superpositions of partial polaritonic waves, 

E(z,t) = ^AjE^t) 

3 

_ 'y ' e i(Kj+nK)z—i(u+nUt) (29) 

nj 

P(z,t) = J2A j P nj e i{k i +nK)z - i{u+nUt \ (30) 

nj 

where E n j and P n j are normalized eigenvectors and Kj 
are eigenvalues described and calculated in Sec. II V Al In 
the vacuum region (z < 0), the electric field is given by 

E(z, t) = e iqoZ e- iut + r n e- iq " z e~ l{uJ+nn)t , (31) 

n 

where q n = {ui + nCl)/c is the wave vector and r n is 
the amplitude of the outgoing Bragg replica n. The am- 
plitude r n is normalized to the unity amplitude of the 
incoming light wave, which is described by the first term 
in the right hand side of Eq. (jolj) . Coefficients Aj, which 
stand for the amplitudes of the partial waves, are found 
from Maxwell's boundary conditions. These require con- 
tinuity of the electric and magnetic fields, E(z,t) and 
H(z,t), across the interface z = 0. Applying the bound- 
ary conditions to the electric field Eq. ([29]) and (l3"Tj) and 
to the magnetic field calculated via dH/dt — cdE/dz 
results in the following set of equations 

g — iut _|_ r ^ e -i(ui+nfl)t _ AjE n je~ i i w + n ^) t 

n nj 
e ~iojt _ ^""^ r ^ e -i(u>+nW)t 

n 

= Y^AjE nj ~ Kj+nK e-^ + ^\ 

Qn 

nj 

which have to be satisfied at any time t. Thus we equate 
coefficients of the same exponentials, which then leads to 
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a set of linear algebraic equations for r n and Aj : 

0n,0 + r n = J2^E nj , (32) 

i 

Sn,o-r n = Y, A i E *i- • ( 33 ) 



Equations (1521 and (|3"3"f are solved numerically, using 
a truncation -rig n ^ tib with the total number of 
Bragg replicas Nb = 2ns + 1. Normally, the number of 
Bragg replicas Nb required for an accurate calculation 
of the reflectivity spectra is much smaller than the ma- 
trix dimension M used in the calculation of the polariton 
quasi- momentum x(iS), Nb <C M. For example, to reach 
the accuracy 10 -7 in calculation of CuCl reflectivity us- 
ing the AW parameters a = 2 meV and v = 100 MHz, one 
needs to take only Nb — 21 and M = 81 into account. 



C. Reflectivity spectra 

In our geometry, the light and the AW propagate 
collinearly and thus there is no way to have any spa- 
tial selection of the Bragg replicas. However, each Bragg 
replica has it own unique temporal evolution as can be 
seen in Eq. (pITj) . Thus, using a frequency selective tech- 
nique one can measure the intensity of the n-th Bragg 
replica in the reflected light, 



Rn 



(34) 



The n-th replica carries the frequency uj + nfl, it is down- 
converted (up-converted) with respect to the frequency 
lj of the incoming light, for negative (positive) n. 

Without such frequency selection, the total but still 
spectrally resolved Bragg reflection can be measured. 
Due to the different time evolution of different Bragg 
signals, there can be two opposite ways to measure the 
reflectivity. One is fully incoherent, in this case the tem- 
poral resolution r of the spectrometer is insufficient to 
see the interference fringes between different Bragg repli- 
cas. In other words, the interference terms fully vanish in 
the measured electric field. This happens when fir 3> 1, 
so that r falls into a nanosecond range. In this case, the 
total Bragg reflectivity is given by 



R 



E 



(35) 



Alternatively, the spectrometer can have enough tempo- 
ral resolution and then the opposite limit fir <C 1 is re- 
alized making the total Bragg reflectivity fully coherent: 



R = 



(36) 



Such regime is achievable when the spectrometer resolu- 
tion is in the sub-nanosecond or even picosecond range. 



We show in Fig. 5 both the coherent and incoherent 
signals, R(uj) and R(u>), in comparison with the zeroth 
replica Rq(uj). Note that Rq is the only component which 
survives in the absence of any acoustic pumping; for 
a = it is shown in Fig. 5 by a red (gray) curve, the same 
in all panels. All other spectra shown in Fig. 5 by black 
curves, are calculated for a = 1 meV. They reveal drastic 
changes compared to the bare (tr = 0) spectrum. By in- 
creasing the acoustic frequency from 50 MHz to 1 GHz, a 
strong spike in the reflectivity R, which is associated with 
the acoustically modified LP branch moves towards the 
Restrahlen band, gets thicker, and finally saturates into 
a cusp-like feature close to the bare TO-phonon transi- 
tion energy. Well below the TO-phonon energy, the spike 
is exclusively due to the first-order acoustically induced 
band gap discussed in Sec. lIII Bl Its thickness is given ei- 
ther by the band gap width Ai [see Eq. (|2"!Zj) ] or by the 
phonon damping, whichever is greater. A similar spike 
which comes from the UP branch is well resolved above 
the Restrahlen band, for the acoustic frequencies between 
100 and 300 MHz. The coherent Bragg reflectivity R has 
similar properties, however the spectral changes are en- 
hanced considerably by pronounced dispersive features 
which appear due to the strong interference effect. 

At higher acoustic frequencies, the contribution of the 
central Bragg replica, R$, becomes dominant and all 
other replicas vanish as v — > oo, so that all three spectra 
R(u>), Ro(u), and R(cv) finally coincide. This happens 
because the light wavelength becomes large compared to 
the acoustic wavelength and thus the AW modulation of 
the semiconductor structure is not resolved by the light 
field. The acoustically modulated semiconductor works 
in this case as an effective homogeneous medium with 
a modified w-dependent dielectric constant. The cusp- 
like structure in the reflectivity which evolves from the 
dispersive (N-shaped) feature in R [see Fig. 5(b)] can 
be understood as a manifestation of the band gaps ac- 
cumulating around the TO-phonon frequency. Indeed, 
according to Eqs. (|20|) and (|2"Tj) . at large values of the 
acoustic wave vector K, there must be a plenty of band 
gaps stacked together around Wt- 

In this high-frequency limit, there are no qualitative 
changes in the reflectivity when the coupling to the AW 
increases, as it is clear from Fig. 6(b), where the reflectiv- 
ity R (R s» i?o) is shown for different values of the cou- 
pling strength a. Indeed, in this case only the strength 
of the cusp, its splitting from the Restrahlen band, and 
the transparency of the latter grow with a. In the op- 
posite case of small acoustic frequencies, the increasing 
coupling strength does lead to qualitatively new features 
in the spectrum. In particular, new peaks correspond- 
ing to the second and third acoustically induced band 
gaps, both in the UP and in the LP branches, grow up 
in v = 50 MHz spectra as shown in Fig. 6(a). 

The influence of the first- and higher-order band gaps 
on different Bragg replicas is analyzed in more detail 
in Fig. 7. The top panels display the polariton quasi- 
momentum dispersion, similar to those shown in Figs. 3 
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FIG. 5: (a) Incoherent and (c) coherent total Bragg reflectivity calculated for a — 1 meV and different acoustic frequencies, 
without and with interference terms, respectively, (b) The same frequency central Bragg replica. The bare reflection spectrum 
of damped TO-phonons (with no AW) is shown by red (gray) thin curves. 



and 4. The Bragg condition is illustrated in Fig. 7 by 
gray horizontal lines (separated by K/2) and by arrows 
all having the same length equal to K. The main spikes 
in the full reflectivity spectrum which are shown 

in the bottom plots, appear exactly at the energies of 
the first-order band gaps. These spikes are due to one- 
phonon transitions between positive->r and negative-^ 
parts of the polariton dispersion h(cj). This effect, which 
is already discussed above in terms of the bare polari- 
ton dispersion uj (k) (see Figs. 1 and 4), is now demon- 
strated taking into account the realistic acoustically mod- 
ified dispersion x{uj) . With increasing acoustic frequency 
Q, and, consequently, increasing acoustic wave vector K, 
the same Bragg condition can only be satisfied at a higher 
frequency. Thus, the band gaps and the spikes are blue- 
shifted, at the same time getting wider due to a larger 



TO-phonon component. This is visualized by two and 
tree times longer arrows in Figs. 7(b) and (c). Below the 
Rcstrahlcn band, these arrows point to one-phonon tran- 
sitions that occur, respectively, at nearly two and three 
times larger frequency. The band gaps are more pro- 
nounced at larger Q and manifest themselves as longer 
and flatter plateaus in the dispersion [see Fig. 7(c), top 
panel]. From the practical viewpoint, the effects demon- 
strated in Figs. 5-7 are nothing else than a very efficient 
way of changing the terahertz properties of semiconduc- 
tor structures simply by tuning the frequency and the 
intensity of the pumping AW. 

The A^-phonon transitions (with N — 2, 3, and 4) are 
shown in Fig. 7 by a series of N arrows of different color. 
They are responsible for higher-order acoustically in- 
duced band gaps and for new spikes, which appear in the 
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FIG. 6: Total Bragg reflectivity R(cj) calculated for the acoustic frequency v — Q/(2n) — 50 MHz (a) and 1 GHz (b), and for 
different values of the coupling strength a (acoustic intensity 7 ac )- Small vertical arrows in panel (a) point at the spike growing 
due the third-order band gap. The bare reflection spectrum of damped TO-phonons is shown by red (gray) thin curves. 



reflectivity when a increases (see also Fig. 6). In particu- 
lar, Fig. 7(a) demonstrates that the spikes due to second-, 
third-, and fourth-order gaps show up in the Bragg repli- 
cas n = —2, —3, and —4, respectively. The higher-order 
Bragg replicas, however, accumulate the contribution of 
all lower-iV phonon transitions in higher orders. The 
higher-order effects of lower- N gaps are sometimes even 
stronger than the contribution of the same AT = \n\ gap 
to the Bragg replica r„, see e.g. r_3 and r_4 replicas in 
Fig. 7 at the frequencies cj « u>\ and 0J2 ■ This effect will 
be discussed in more detail in Sec. lIVDl 

The incoming light beam excites polaritons with posi- 
tive values of the quasi-momentum k > 0, and then these 
polaritons are scattered by the AW to the negative- x part 
of the polariton dispersion. Such scattering is accompa- 



nied by emission of acoustic phonons and consequently 
by down-conversion in frequency. These down-converted 
polariton waves correspond to negative-n components of 
the electric field, both in semiconductor and in vacuum, 
see Eqs. (|29[) and (|3"Tj) . That is why all spikes in the 
reflection corresponding to the AW-induced band gaps, 
appear only in Bragg replicas having negative n and no 
such features can be seen when n is positive. 

By increasing the coupling strength a, A^-phonon tran- 
sitions responsible for the band gaps reflectivity spikes 
become increasingly dressed by higher order transitions. 
In particular, the A^-th order gap is formed due to multi- 
ple processes of simultaneous emission of N + s phonons 
and absorption of s phonons, where s can be any integer. 
In the present calculation however, with the parameters 
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FIG. 7: Acoustically modified polariton dispersion of damped TO-phonons (top panels), contributions of different Bragg replicas 
\r n \ 2 with n = 2, 1, 0, -1, -2, -3, and -4, as well as the total reflectivity spectrum R(uj) = ^2 n \r n \ 2 (bottom panels), calculated 
for a — 2meV and v — Q,/(2n) — 50 MHz (a), 100 MHz (b), and 150 MHz (c). Arrows and vertical lines indicate positions of 
the band gaps due to iV-phonon dressed transitions. The acoustic wave vector takes the values of (a) K = Kq, (b) K — 3Ko/2, 
and (c) K = 2Kq, where Kq = 0.1555 /Urxi - . 



12 



used throughout the paper, this number never exceeds in 
reality |s| = 20, even for very large acoustic intensities 
(up to a = 3 meV) . 



D. Bragg replicas: Analytic approach 

At low acoustic intensities (small a), the Bragg replicas 
can be calculated using the perturbative approach intro- 
duced in Sec. lIIIBl In this approach, the A-th pertur- 
bation order describes the acousto-optic effect with par- 
ticipation of A phonons (no higher-order dressing is in- 
cluded) , and in particular the opening of the gaps A n due 
to A-phonon transitions. Below we calculate acousto- 
optic effect in each A-th Bragg replica in its lowest (A- 
th) order, calculating both up-converted r + jv and down- 
converted r_jv replicas, and show explicitly that the ef- 
fect is resonantly enhanced only in the down-converted 
replicas r_i, r_2, etc. 

The perturbation-theory result for the A-th Bragg 
replicas has the following recursive form (for derivation 
see Appendix |D|) : 



r±N 



s=l 
N 



a ±n = ~y^ J A ±(N _ s) E± s 1 + ^ =s , (38) 
N 1 

E ±N = a^^n 



(39) 



a±N fj[ D± s 

where D n (uj, k) are defined in Eq. ([T5j) . 

fk + nK\ 2 „ k + nK 
a n = - 1 , p n = , (40) 

\ p J q 

(we have neglected here any dependence on ft), q = uj/c, 



and p = Jelq. The starting A = values in the recur- 
sion are given by 



''o 



l-A) 
l + /3o 



A 



1 + A> 



E 



1 



(41) 



To a first approximation, the quasi-momentum k in the 
above equations is taken in zeroth order, i.e. is calculated 
according to the bare polariton dispersion 



jfe°(w) 



(42) 



which satisfies Eq. (fT!?]) and is inverse to w = u>°(k). Us- 
ing this bare dispersion relationship, Do(u>,k) = 0, the 
functions D n (uj,k) are simplified to 



a n D n = -2. 

2uj t 



lo% (k + nK) 2 



k 2 



P 



(43) 



where any dependence on Q is neglected. 

The explicit form of Eq. (|43|) allows us to analyze the 
resonances seen in Bragg replicas. Since D n stands in the 
denominator of the electric field components Eq. (f39| , r n 
have maxima when a n D n is minimized. This happens 
when (k + nK) 2 ss k 2 , leading to a Lorentzian spectral 
line shape peaked at the Bragg condition, 



„ nK 
Re(fe) = -— ; 



(44) 



and having the widths proportional to Im(k) which at 
this point is fully determined by the TO-phonon damp- 
ing 7. The incoming light excites polaritons with positive 
wave vectors Ke(k) > 0, thus the Bragg condition is ful- 
filled only for negative n < 0, so that the peaks in the 
reflectivity are observed in the negative (down-converted) 
Bragg replicas only. 

The amplitudes of the three lowest-order Bragg repli- 
cas have the following explicit form 



K 



r-i 



r-2 



r-3 



on 



2an 



q (1+A,) 2 ' 
, K 2 

'7 (i + A)) 2 

K 2 



1 + /3- 



a 3 - 



q (1 + A)) 2 



3a 3 + 2/3_!+^_ 2 



a_ 3 D_ 3 Z)_2i?-i 



1 + Po 



1 + /3-1 
l + /3o 



a_ 



"0 
,D 3 



(45) 
(46) 

(47) 



The analytic result Eqs. P5 |) - (j4T|) is demonstrated in 
Fig. 8(a) in comparison with the full numerical calcu- 
lation; both are done for a moderate coupling strength 
a = 0.4 meV. The agreement is good everywhere except 
the spikes, in which only positions are reproduced well, 



see a zoomed area in the inset of Fig. 8(a). Indeed, with 
purely Lorentzian peaks in the spectra, the magnitudes of 
the replicas sometimes exceed unity which results in fully 
unphysical reflectivity. This artefact is however removed 
below and the spectrum is refined by taking into account 
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FIG. 8: Bragg replicas \r n \ calculated numerically (black 
lines) and analytically [red, green, blue, and magenta (gray) 
lines] for (a) n = —1, —2, —3, and (b) n = 1, 2, 3, and 0. 
Parameters used in both calculations: the coupling strength 
a = 0.4 meV, the acoustic frequency v = 50 MHz, the TO- 
phonon damping 7 = 0.2 meV. The inset zooms in on the 
spectral region of the first acoustically-induced band gap and 
spikes in the Bragg reflection due to one-phonon transitions. 



the acoustically modified polariton dispersion. For pos- 
itive values of n, the replicas are shown in Fig. 8(b) in 
order to demonstrate that no spikes at the energies of 
AW-induced band gaps are seen either in analytical or in 
numerical calculation. 

The first down-converted Bragg replica r_i has a spike 
at the frequency oj± of the lowest acoustically-induced 
energy gap, see Eq. (f^Tj) . The spike corresponds to a 
one-phonon transition between polariton states in the 
positive-A: and negative-fc branches of the dispersion, as 
shown in Fig. 7 and discussed earlier in Sec. lIV CI Note 
that although we concentrate here mainly on the LP part 
of the spectrum, similar spikes can be seen also in the 
spectral region of the UP branch. 

In the second replica r_2, resonances appear at the 
energies of both first and second band gaps, uji and u>2- 
The first term in square brackets in Eq. (|4"f5|) describes a 
composite two-phonon process with a virtual intermedi- 
ate state having only one phonon emitted. That is why 
the first term in Eq. (|4"fJ)) shown in Fig. 9 by a dashed 
red line is peaked at both energies wi and u>2. The fi- 
nal state has the dominant role and so the resonance at 
uj ~ L02 is much stronger. The second term shown by blue 
(gray) line in Fig. 9 has only one peak at U) rs u>i (in the 
LP part of the spectrum) as it describes a one-phonon 
contribution in second order. This one-phonon process 
is the only component of the first-order Bragg replica 
r_i. Surprisingly, its dominant role is well seen in r_2 
as well, see the sum of both terms, the black solid line 
in Fig. 9. This physical picture helps us to understand 
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similar features in higher Bragg replicas. In particular, 
the first term in Eq. (|47[) is due to a three-phonon tran- 
sition having as intermediate (virtual) states those with 
one and two phonons emitted. The last term is the one- 
phonon process taken into account in third order, and 
the middle term is a mixed contribution of both one and 
two phonon-assisted transitions. Consequently, the reso- 
nances in r_3 occur at u> « oj±, 0J2, and W3 and are due 
to different processes with participation of three acoustic 
phonons. 

Let us now refine the analytic calculation in order to 
get rid of the artificially high analytic reflectivity which 
is detailed in the inset to Fig. 8(a). To do this, we sim- 
ply replace the bare polariton dispersion k = k (ui) in 
Eqs. (l40l) and (|43| by the acoustically modified dispersion 
k = k(lo) calculated in Sec. lIV Al [see Eq. ([28]) ]. Conse- 
quently, a modified Bragg condition which follows from 
Eq. (|43[) and determines the acousto-optic resonant ener- 
gies ui n becomes 



Re[x(w|„|)] 



nK 



(48) 



[cf. with Eqs. (PHI) and (|2T|) ]. The refined Bragg condition 
Eq. (|48)) has been used already in our analysis of the spec- 
tra shown in Fig. 7. From the viewpoint of the perturba- 
tion theory developed in the present section, the calcula- 
tion of the modified dispersion x(oj) can intrinsically take 
into account a large number of higher perturbation or- 
ders. Thus, it is not fully consistent with the calculation 
of r n limited to the lowest order |ro| only. From the prac- 
tical viewpoint, however, such a refinement works per- 
fectly well for sufficiently narrow resonances, as can be 
seen in Fig. 10(a) in which the spectral lines are shown by 
the solid red and green (gray) curves in comparison with 
full calculation (black curves). Both the height and the 
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FIG. 10: Bragg replicas |t"_i| 2 and |r*_2 1 2 calculated ana- 
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eters are the same as in Fig. 8. The double-side arrow shows 
the width of the first acoustically induced gap calculated via 
Eq. (H2J) 



V. DISCUSSION 

A. Low and high acoustic frequency limits 

The terahertz properties of semiconductors and their 
changes due to varying acoustic frequency f2 has been al- 
ready discussed in Sec. lIV Cl In particular, Fig. 5(a) has 
demonstrated an efficient manipulation of the spectral 
properties, changing both the positions of the spikes and 
their widths simply by tuning tt. Moreover, in the lim- 
iting case of very large f2, the wavelength of the spa- 
tial modulation caused by the AW turns out to be much 
smaller than the wavelength of the propagating light, so 
that the modulation does not produce any diffraction. As 
a result, all acoustically induced spectral changes are con- 
centrated in the central Bragg replica Rq [cf. Figs. 5(a) 
and (b)], while the contribution of all higher replicas van- 
ishes at — > oo. This is more clearly demonstrated in 
Fig. 11 where the lower panels shows diminishing contri- 
bution of the higher (\n\ > 0) Bragg replicas to the total 
reflection. The results are summarized in Fig. 12 where 



width of the spikes are reproduced well. The narrow and 
unrealistically high Lorentzian lines shown by the dashed 
curves are fully determined by the damping 7. In the re- 
fined analytic calculation, the lincwidth is determined by 
the band gap width, which is properly incorporated into 
the modified dispersion yt(uS). A slight difference between 
the numerical and refined analytical calculations is due 
to the fact that for such a moderate acoustic intensity 
(a = 0.4 meV), higher-order corrections (with s > |n|) 
to the Bragg replicas r„ which not included in Eqs. ([571) - 
are already perceptible. They become crucial in case 
of much wider gaps in the UP part of the spectrum, see 
Fig. 10(b), where the improvement due to the dispersion 
refinement is rather minor. 

Finally, a calculation of the n = Bragg replica, which 
has been done up to now only in zeroth order, can be fur- 
ther improved by inclusion of second-order corrections: 



ro = 



I-/80 
1 + A> 



(1 + A)) 3 ai£>ia_i£>_i 



(49) 



Clearly, this expression contains resonant features due 
to the D-i in the second term. However, unlike the 
higher order replicas which exhibit the Lorentzian line 
shape, this resonant term is added to the background 



reflection, Eq. (|41j) . Thus, \tq\ has instead a more pro- 
nounced dispersive behavior, due to the mixing of the 
first and second terms in Eq. (l4*9"l) . Such dispersive fea- 
tures manifest themselves as N-shaped resonances clearly 
seen in Fig. 5(b) and become more and more prominent 
with increasing acoustic frequency (see also n = repli- 
cas in Fig. 7). 
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FIG. 11: The total contribution of higher Bragg replicas, 
X] n ^o l r ™| 2 ' f° r t ne coupling strength a = 1 meV and different 
acoustic frequencies. 
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FIG. 12: The spectral area of higher Bragg replicas, calcu- 
lated according to Eq. (|50[) for a = 1 meV and different values 
of the acoustic frequency v — Q/(2n). 
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the spectral area of the higher replicas is plotted against 
f2. This spectral area is calculated according to 



^2 KM| 2 du . 



(50) 



Surprisingly, a similar effect takes place in the opposite 
limit of low acoustic frequencies. Although the reduction 
of the contribution of higher replicas in this case is not 
very well illustrated by Fig. 11, and technically this limit 
is hard to achieve with the present numerical method, 
an overall trend is such that the spectral area is get- 
ting smaller and possibly vanishes at fl — > 0, as seen in 
Fig. 12. In this case, the short-period acoustic modula- 
tions work again as a homogeneous dielectric medium, 
which can be be characterized by an effective (modified) 
dielectric constant. In a sense, this picture is similar to 
the effective mass approximation for electrons, valid in 
bulk semiconductors as well as in short-period semicon- 
ductor superlattices where the Bragg scattering can also 
be neglected. 

A clear bump just above the Restrahlen band seen in 
the top spectrum in Fig. 11 which is calculated for a rel- 
atively low acoustic frequency is due to accumulation of 
band gaps near the bottom of the UP branch (a similar 
accumulation of gaps takes place at — > oo near the 
top of the LP branch; it has been already discussed in 
Sec. lIV Cl) . This narrow spectral region remains infrared- 
active in higher Bragg orders due to existence of polari- 
tons which have momenta the same order of magnitude 
as the acoustic wave vector K (very small in this case) 
and thus can be scattered by such an AW. 



B. Electric field distribution and interaction length 

When the wavelength of the incident light and the pe- 
riod of the acoustic modulation are comparable as shown 



FIG. 13: The electric field intensity \E(z,uj)\ in vacuum 
and semiconductor regions, calculated for a — 1 meV, v = 
50 MHz, and different frequencies of the incoming terahertz 
light. Black in the logarithmic gray scale map corresponds to 
|£| 2 ^0. 



in Fig. 13, it results in significant spectral changes which 
are mainly due to a strong Bragg scattering of the po- 
lariton waves from the propagating periodic patterns in- 
duced by the AW. The electric field distribution in Fig. 13 
corresponds to the terahertz spectrum shown in the top 
panel of Fig. 5(a). While in the vacuum region the modu- 
lations of the electric field are determined exclusively by 
the light wavelength (resp. lj) and are caused by a su- 
perposition of the incoming and outgoing waves. In the 
semiconductor region the period of modulations is the 
same for all uj and is given by the acoustic wavelength 
2n/K. At around Tho — 5meV the AW produces a half-A 
modulation (A is the wavelength of the light propagating 
inside semiconductor) that results in a destructive inter- 
ference of co- and contra-propagating polariton (light) 
waves, and thus in quick attenuation of the electric field 
seen as a narrow black stripe in Fig. 13. 

The decay length of the electric field gives us an upper 
bound for an estimate of the acousto-optic interaction 
length Z; nt . In resonant acousto-optics, such a length 
characterizes the scale in real space where the light ef- 
fectively couples to the ultrasonic wave by means of the 
TO-phonons. The interaction length calculated at the 
resonant frequency u = uj\ of the first AW-induced band 
gap is evaluated analytically in Appendix [El assuming 
-K7int 3> 1- Taking into account the finite TO-phonon 
damping 7, it takes the following form 



In 



Kc 2 (ujf - uj 2 ) 2 + 4 7 2 w 2 



(51) 



ae b ix> t 

Figure 14 shows the interaction length numerically eval- 
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FIG. 14: Interaction length, corresponding to the first 
acoustically- induced LP band gap u) = OJi, calculated analyt- 
ically via Eq. (|51[l and evaluated numerically from the electric 
field profile, for a = 1 meV and different values of the acoustic 
frequency (top axis) or the detuning u>t — wi (bottom axis). 
The dashed line shows the acoustic wavelength. 



uated from the decay of the electric field at w = w\, 
in the LP part of the spectrum. It also shows an ex- 
cellent agreement with the analytic result Eq. (|51[) , over 
the whole range of acoustic frequencies ft, or detuning 
5u> = u>t — uj\ . When decreasing the detuning, the inter- 
action length also decreases monotonously and at very 
small values of Suj is limited by the TO-phonon damping. 
It exceeds the acoustic wavelength 2ir/K (also shown in 
Fig. 14 by a dashed line) by a few tens to hundred times 
only. 



VI. CONCLUSIONS 

We have studied, both numerically and analytically, 
the far-infrared properties of acoustically pumped semi- 
conductor crystals, in a geometry of collinear propaga- 
tion of the light and ultrasound waves. Owing to the 
cubic phonon anharmonicity, TO phonons mediate the 
coupling between the light and acoustic fields that results 
in drastic changes of the terahertz properties of semicon- 
ductors and in enhancement of the acousto-optic effect 
by several orders of magnitude compared to the conven- 
tional acousto-optics. 

We have developed an efficient numerical method for 
calculation of acoustically induced band gaps, electric 
field distribution, and such infrared properties of semi- 
conductors as transmission, reflection, and Bragg scat- 
tering of light. The method is based on the plane waves 
expansion of the electric field and TO phonon-induced 
polarization. It allows us to solved the coupled Maxwell 
and material equations (including a phonon anharmonic 
source term due to a pumping AW) by mapping them 
into a linear matrix problem. We have provided a de- 



tailed analysis of the Bragg replicas of the incoming 
light (the Bragg replicas are frequency-selected in the 
present case of collinear geometry) and demonstrated a 
one-to-one correspondence between the spikes in the re- 
flectivity and the acoustically induced gaps in the polari- 
ton quasi-energy spectrum. The gaps open at the ener- 
gies u>m satisfying the Bragg condition Eq. (l4"51) . manifest 
themselves as spikes in down-converted Bragg replicas r n 
with n = — 1, —2, . . . , — N, and correspond to 7V-phonon 
transitions (dressed with higher orders) between different 
polariton states in the crystal. In particular, the reflec- 
tivity is enhanced considerably when the acoustic wave- 
length is close to half integer (N/2) wavelength of the 
light in crystal. In this case a polariton wave is scattered 
by the periodic pattern of acoustic modulation of the 
TO-phonon field emitting A?" acoustic phonons. The co- 
propagating polariton wave and contra-propagating scat- 
tered wave interfere destructively to produce the elec- 
tric field decay into the crystal. We have calculated the 
length of such decay which, in turn, limits the inter- 
action length between the light and ultrasound waves. 
The latter can be as short as only tens of the acoustic 
wavelengths, much shorter than in conventional acousto- 
optics. 

Our analytic approach is based on the perturbation 
theory developed in the limit of weak coupling a ju t <C 1 
and allows us to evaluate and analyze the band gap 
widths, the Bragg replicas, and the interaction length. 
Moreover, being in good quantitative agreement with nu- 
merical results, the analytic approach remains a reliable 
and accessible tool for calculation of the terahertz prop- 
erties up to a moderate coupling strength (a = 1 meV 
in CuCl, corresponding to the acoustic intensity / ac = 
14.8kW/cm 2 ). We have also provided a refinement of 
our analytical method which uses an acoustically modi- 
fied polariton dispersion and considerably improves the 
calculation of the sharp spikes in the reflectivity. 

Finally, in the limits of very low and very high acoustic 
frequency ft, the contribution of all higher Bragg repli- 
cas r n with n ^ vanishes, but the central replica ro still 
shows considerable spectral changes which grow dramat- 
ically with the acoustic intensity (or with a). In these 
limiting cases, the wavelength of the propagating electro- 
magnetics waves is either too small (ft — > 0) or too large 
(it — >• oo) compared to the acoustic wavelength to feel the 
periodical modulation of the crystal, therefore, no Bragg 
scattering occurs. Instead, the polariton waves propagate 
like in a homogeneous medium which can be character- 
ized by an w-dispersive effective dielectric constant. 
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Appendix A: Fourier transform of the electric field 
and polarization 

In this section, we show that the partial waves of 
the electric field and TO-phonon polarization, which are 
found as eigenfunctions of Eqs. (|5|) and (j6|), can be writ- 
ten in the form of a single Fourier transform, Eqs. (JjfJ and 
©. Due to the periodicity introduced by the AW, each 
partial wave can be expanded into Fourier series, both in 
time and space. In problems similar to the quasi-energy 
spectrum calculation where the wave vector k is a fixed 
real- valued parameter, it is convenient to use the follow- 
ing space-time Fourier transform of a partial wave: 

E(z,t) = J2£ nm e iik+nK)z e- i{u+nn+mn)t , (Al) 



pom) = J2 V " 



(A2) 



Substituting the above expansions into Eqs. J3]) and © 
and equating coefficients at the same exponentials results 
in the following eigenvalue problem 



(Jfe + nKf - %(u) + nil + mil) 2 

An 



n 



-y (w + nil + mil) z V nm , (A3) 



ui 2 — 2ij(oj + nil + mil) — (cj + nil + mil)' 



s tu 2 

47T 



(A4) 



which has to be solved to find eigenvalues U) and eigen- 
vectors £ nm (and V nm ). It is convenient however to in- 
troducing a new eigenvalue <D = cj + mil which makes the 
eigenvalue problem to- independent, so that the eigenvec- 
tors £ nm an d £nm' calculated for different to and to' turn 
out to be proportional to each other and corresponding to 
the same, TO-independent eigenvalue Co. In other words, 
all terms in Eqs. (IA1|) and (|A2[) having the same index 
n are linearly dependent, which effectively eliminates the 
TO-summation bringing the above series to the form of 
Eqs.© and ©. 

Note that the second index to of the formal double 
Fourier transform £ nm and V nm has nothing to do with 
the index j in E n j and P n j used throughout the pa- 
per to label different partial waves (normalized polari- 
ton eigenstates) while the first index n denote different 
Fourier/Bragg components of the same partial wave. 



(i) Let us first consider the quasi-energy dispersion 
u)(k) in which a real- valued wave vector k is a given num- 
ber (parameter). We assume 7 = as done in Sec. lIII A"l 
To linearize the eigenvalue problem, we introduce two 
new vectors, 

A n = (io + il)P n , (Bl) 
B„ = {uj + il){e b E n +A^P n ) 

= {u + tt)e b E n + A-nA n , (B2) 

and substitute them into Eqs. (JSJ) and (fTTj|) . The eigen- 
value problem then takes the form 



uP n = -nilP n + A n , 

TP OP _L Bn 4lrAn 

ujE n = —n\lhj n -\ , 

LuB n = -~nilB n + c 2 (k + nK) 2 E n , 



(B3) 
(B4) 
(B5) 



UjA n 



-nilA r , 



47T 



En 



+u 2 P n + 2aw t {P n+1 + P„_i) . (B6) 



The above set of equations can be written in a compact 
matrix form: 



wX = WX , 



(B7) 



where W and X are, respectively, AM x AM matrix and 
AM- vector defined as (the truncation of the matrices the 
and the definition of M are introduced in Sec. lIII Al 



-iln 


i 



0J 2 1 + 2<TLJ t f 



-iln 






-iln 






c 2 (fci + Kn) 2 


—iln 



X 



(B8) 



(B9) 



1, h, and / are M x M matrices having the following 
properties: 



(nP)„ = nP n , 

(fP) n = P n ^+P n+1 , 



(BIO) 
(Bll) 
(B12) 



Appendix B: Linearization of the eigenvalue problem 



The eigenvalue problem given by Eqs. © and ([TU]) is 
nonlinear, both in uj and k, and is solved either to find 
the polariton eigenfrequencies in Sec. lIII Al or to find the 
complex polariton wave vectors in Sec. lIV Al Let us show 
that in both cases this set of equations can be reduced 
to a standard linear eigenvalue problem. 



and A, B, P, and E are M-vectors. 

(ii) In the opposite case of the eigenvalue problem to 
be solved for the complex- valued momentum k, keeping 
the real frequency to constant, we do essentially the same. 
However, before introducing a new vector, we express the 
polarization P n from Eq. (|10[) which can be written as 



E G 



P — E 

run 1 m ' -^n j 



(B13) 
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where 



4-7T 



lu 2 — 2i~f(uj + nft) — (lu + nfi,y 



+2<TLUt(S n ,m+l + 5 n ,m-l) 



(B14) 



Noting that G nm does not depend on k (it depends on 
u only), this matrix can be inverted that results in the 
following equation 



4?r 



(u + nn fP n = HnmK 



where the matrix H has the form 



47T 



H nm = —(u + nQ) 2 (G~ 1 ) r 



(B15) 



(B16) 



Comparing Eqs. © and (|B15|) and also introducing a new 
vector F defined as 



, u + nf2 s , 

F n =(k + nK + ^/T b ) E n , 



(BIT) 

a linearized matrix eigenvalue problem takes the form 

VY = fcY, (B18) 

where V and Y are, respectively, 2M x 2M matrix and 
2M-vector, defined as 



-nK 



£6 ' 

H 



ui + hQ 



1 



-hK + ^/efc 



wl + fiQ, 



Y 



(B19) 
(B20) 



The procedure described in Sec. lIV Al to refine the po- 
lariton wave vectors kj^to using the best calculated cen- 
tral root >c(uj), utilizes the translational property of the 
modified polariton dispersion Eq. (f24)) and the following 
linear interpolation between neighboring points in fre- 
quency, U>i ^ LU — jCl ^ LOi+1'- 



kj(u>) = jK + x(Lo—jQ) 



(B21) 



x(u i+1 ) - x(c^) 
= a{u)i)-\ (w-jSZ-w,). 

Appendix C: Analytic evaluation of the band gaps 

When calculating the band gap widths Ajv we have to 
take into account the effect of mixing (and anticrossing) 
of two polariton waves. One wave has the wave vector 
k = kf4 = NK/2 and refers to the n = replica, the 
other has the wave vector k — k^ — NK = —NK/2 and 



refers to n — —N component of the electric field. To take 
them into account in the lowest (iV-th) order we truncate 
the tridiagonal matrix equation Eq. (fT4|) to 



= 0, (CI) 



Do -a 


. 


. 




a E 




-a D-i 


—a . 


. 






l 


-a 


D-2 ■ 


. 




a- 2 E- 


2 





. 


■ D_ N 






JV 



where D n = D n (kN,u) and 



D n (k,u) = — 

ZUt 



2 2 
LU, — LU — 



(k/q) 2 1) 



(C2) 

BbLu/c. 



neglecting 7 and O [cf. with Eq. ([To])], and q 
We then solve Eq. (|C1I) and find small deviations of u 
from the bare gap positions u^ given by the equation 
Do(k]y ,lun) — 0. First of all, we find deviations of Do 
from zero by solving 



det 



D -a 
— a D-i —a 
-a D-2 







D, 



-a L>_i 







-a D 



= 0. 



(C3) 



where the symmetry property D s = Dn- s is used. This 
results in the following quadratic equation for Dq: 



DZSu - 2D <T 2 S 12 + (t 4 5 2 2 = . 



(C4) 



where (for the fixed N) 



D- n -a 
-a £>_ n _i 



Sn 



det 





— a 

D- n -2 



D- m -i —a 
-a D- m 



Noting that 

S12 — S11S22 = C (S23 — S22S33) = • ■ ■ 

we find 

d o[Un) = c 

On 



= (j2(JV-2) 



(C5) 
(C6) 

(C7) 



where the two values are the energy positions of the 
top and the bottom of the iV-th gap. The gap width 



LU 



N 



D (lu n ) - D (u N ) = 



u N \ is then proportional to 

2a N 2a N 



Sn 



D\D 2 _ 2 ...D S 



(C8) 
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where we have neglected all orders higher than N. In the 
last equation, s — 1 and [^] = -y for even TV, and s = 2 
and [-y] = for odd AT. Finally, the band gap A^v is 
evaluated using the linear approximation 

dD 



D (u+) - D (u N ) « (ut - u N ) 



duj 



(C9) 



where the derivative is calculated and simplihcd to 



dD _ w|tj| + (tj| - w 2 ) 2 

diO L0 t LOUJ R 



(CIO) 



using the condition D${k n ,uj) = for uj = lun, taken on 
the bare polariton dispersion. 



Appendix D: Bragg replicas in the limit of small 
perturbations 



Substitution of the plane wave expansion Eqs. (J29J) and 
(|30|) into the Maxwell and material Equation ([5]) and ([6]) 
results in the following tridiagonal matrix equation for 
the electric field components 

& n jD n jE n j — a(a n+ ijE n+ ij + t jE n -\ j^j = , 

(Dl) 

where 



D nj = D n {Rj,u) , 
ot n j = a n {kj,uj) , 



(D2) 
(D3) 



and D n and a n are given by Eqs. ([15]) and (fl6|) . respec- 
tively. Equation (|D1[) is similar to Eq. (TUT) but now in- 
cludes explicitly the index j which refers to a specific 
partial polariton wave. Using the translation property 
Eq. (J2U) we obtain 

D n; j +s (oj) = D n+Si j(u- sQ) , (D4) 
a nJ+s (uj) = a n+Sij (u - sfl) , (D5) 

and then, neglecting the acoustic frequency J7, we are 
able to introduce universal vectors E n , D n , a n , and k n : 



E n +j(oj) — E n j(u>) , 

D n+j (uj) = D nj (uj) , 

a n +j(uj) = a nj (uj) , 

Hj(io) — k{uj) + jK . 



where x(uj) is the central root introduced in Eq. (|28 
Then Eq. ((fTTT) takes the form 



(D6) 
(D7) 
(D8) 
(D9) 



a n D n E n - <7[a n+1 E n+1 + a n -iE n -ij = (DIO) 

which is formally equivalent to Eq. (fn| but has a differ- 
ent meaning. This is because it effectively couples com- 
ponents of different partial waves (having different j). In 



the following we drop all tildes in the above quantities 
for simplicity of notations. 

To solve this equation perturbatively we note that 
E n = 0(a' n '). Neglecting higher orders in a, we obtain 
from Eq. (jDTT)|) 

a\D\Ei = aa^Eo , 
a 2 D 2 E 2 — croi\Ei , 



q±(jv-i) 
a±NE>±N 



or, more explicitly, 



E± N = v 



N a 



N 1 



(Dll) 



(D12) 



where we have used the starting value Eq = 1. 

To calculate the reflection, we use the boundary con- 
ditions Eqs. (|3"2"j) and ([55)1 which take the form (also re- 
placing Aj —> A-j): 



(D13) 
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n,0 



Y,A^E n+j (l + p n+j ), (D14) 



where 



x{u) + nK u 

Pn{u) = 1—, , q{oj) = - . (D15 

q(u)) c 

In the perturbation approach, Aj and r n have the same 
properties as E n , namely Aj — 0{a^) and r n = 0(cr' fl ' ). 
We calculate them in a similar way, keeping only the 
leading-order terms: 



i + /V 

^£0(1 + AO = -AoE^i + px), 

A 2 E (1 + M = -A E 2 (l + p 2 )-A 1 E 1 (l + p 1 ) 



A±n = -^V-^llT 1 ' ( D16 ) 

1 + po 



as follows from Eq. (ID14I) , and 

r = A £/ - 1 = , 
1 + po 

iV JV 

(D17) 

as follows from Eqs. (|D13I) and (|D16I) . 



s=0 
AT 



2^A HN _ s) E ±a ^^- 
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Appendix E: Analytic evaluation of the decay 
(interaction) length 

To calculate the interaction length at the energy of the 
first acoustically-induced band gap w = wi, we take into 
account, in the resonant approximation, the contribution 
of only two co- and contra-propagating polariton waves. 
Keeping in Eq. (|7|) only n — and n — — 1 terms, 



E(z,t) = e 



E e ikz + E_ ie i(k - K)z 



(El) 



we evaluate the imaginary part Q of the polariton wave 
vector 



K 



+ iC 



(E2) 



again, neglecting ft and 7. Then Eq. (fT4|) yields 

chqDoEo — aa-iE-i = , 
a-xD-iE-x — aaoEo = , 



from which we obtain 



D D_x = a 2 



(E3) 



(E4) 



Using the property D_i(fc,cj) = D (k*,uj) = Dq(1c, w) 



valid for the k given by Eq. (|E2j) and introducing the real 



and imaginary parts of Dq as Dq = D' Q + iDq, Eq. (|E4|) 
transforms to 



D$+D>»=<t. 



(E5) 



Since D' = at the gap position u> = u>±, the above 
equation results in 



^'(fc, Wl ) = a 2 , 

or in a more explicit equation for Q: 

,,,2 

Im 



HC/2 + iQ y j 



2LUt<7 . 



(E6) 



(E7) 



where q = y/iEbUj/c. Finally, assuming Q -C if leads to 
_ 2aw t [(if/2 g ) 2 - l] 2 _ 2auj t u J 2 R q 2 



K/q 2 



(E8) 



The interaction length is then evaluated from the decay 
of the electric field intensity and has the form 



hnt " 2Q ■ 



(E9) 
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